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We investigate the influence of quantum fluctuations upon dipolar Bose gases by means of the 
Bogoliubov-de Gennes theory. Thereby, we make use of the local density approximation to evaluate 
the dipolar exchange interaction between the condensate and the excited particles. This allows to 
obtain the Bogoliubov spectrum analytically in the limit of large particle numbers. After discussing 
the condensate depletion and the ground-state energy correction, we derive quantum corrected 
equations of motion for harmonically trapped dipolar Bose gases by using superfluid hydrodynamics. 
These equations are subsequently applied to analyze the equilibrium configuration, the low-lying 
oscillation frequencies, and the time-of-flight dynamics. We find that both atomic magnetic and 
molecular electric dipolar systems offer promising scenarios for detecting beyond mean-field effects. 

PACS numbers: 03.75.Hh,03.75.Kk 



Introduction - Bose-Einstein condensates (BEC) with 
the anisotropic and long-range dipole-dipole interaction 
(DDI) have received much attention, specially after the 
condensation of 52 Cr in 2005 fl[. This investigation pio- 
neered a series of experiments which led to a robust un- 
derstanding of the DDI on a mean-field level. Experimen- 
tal successes include the direct observation of the DDI in 
the time-of-flight (TOF) dynamics [HQ, t ne stabilization 
of a purely dipolar gas [J], and the observation of a d- 
wave Bose-nova explosion [5J. In the meantime, the DDI 
has also been observed even in 87 Rb Q and evidences 
of it have been found in 7 Li Q. Recently, an important 
experiment has been realized, in which the influence of 
the DDI upon the oscillation frequencies of 52 Cr has been 
studied Q. Parallel to these experiments, much theoret- 
ical work has been pursued For instance, building on 
the previous construction of the dipolar pseudo-potential 
[n|, a solution of the mean- field Gross-Pitaevskii (GP) 
equation was obtained 
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[12[, which accounts quanti- 
tatively for the static as well as the dynamic properties 
of trapped dilute 52 Cr gases 0, SB- Nowadays, dipolar 
interactions in magnetic systems are, therefore, consid- 
ered to be relatively well understood in terms of the GP 
mean-field theory. Nonetheless, highly magnetic atoms, 
such as dysprosium [Tj|, or strongly polar heteronuclear 
molecules, exemplified by 40 K 87 Rb |l4-|l6l]. are expected 
to push the understanding of dipolar systems beyond the 
edge of mean-field theory's domain of validity. 

Already in the case of pure contact interaction, ultra- 
cold quantum gases have been the scenario for investigat- 
ing many-body theories. Theoretically, much work has 
been carried out based on the Lee-Huang- Yang (LHY) 
correction for the equation of state 17], when imple- 
mented in terms of superfluid hydrodynamics fl8L [l9j , or 
on quantum Monte Carlo simulations [2(|. Experimen- 
tally, quantum fluctuations have been studied in setups, 



which favor their enhancement, such as in optical lattices 
2l| . or in the presence of Feshbach resonances in both 
bosonic [22| and fermionic [23| systems. In particular, the 
importance of the latter experiment should be stressed, 
as it achieved a precision level which is capable of dis- 
tinguishing the predictions of Bardecn-Coopcr-Schrieffer 
(BCS) mean-field theory from QMC results for the os- 
cillation frequencies of low-lying modes at the BEC side 
of the BEC-BCS crossover. This demonstrates unequiv- 
ocally that collective modes provide an adequate testing 
ground for many-body theories of interacting quantum 
gases. In this letter, we theoretically study quantum 
fluctuations in dipolar Bose gases in the large particle 
number regime. Within the realm of the Bogoliubov-de 
Gennes (BdG) theory, we first calculate the condensate 
depletion and the ground-state energy. Then, by apply- 
ing superfluid hydrodynamics, we discuss the quantum 
corrections to both static and dynamic properties of the 
system focusing on the observability of these effects. 

BdG Theory - Consider a Bose gas of N dipolcs 
with mass M in a cylinder symmetric trapping potential 
f/ tr (x) = Mu? x [x 2 + y 2 + X 2 z 2 )/2 with the trap aspect 
ratio A. The corresponding Hamilton operator reads 



h - 



d 3 



* t (y)^„t(x-y)*(y) 



(1) 

with h = -fr 2 V 2 /2M+ J7 tr (x). Moreover, V&^x) and 
^(x) denote the usual creation and annihilation field 
operators, respectively. Considering the dipoles to be 
aligned along the z axis, the interaction potential reads 



Mnt(x) = 
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where a s = Mg/4irh 2 is the s-wave scattering length and 
add denotes a length scale which characterizes the dipo- 
lar interaction, irrespective whether it is of magnetic or 
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electric nature. At this point, it is useful to introduce the 
dimensionless relative interaction strength edd = Odd/as- 
In order to study dipolar Bose gases beyond the 
mean-field approximation, we follow the Bogoliubov pre- 
scription and decompose the field operator according to 
vE'(x) = ^(x) + Sip (x), where ^(x) is the condensate 
wave function and 5^(x) accounts for the quantum fluc- 
tuations. By considering only the leading term in this 
expansion, one obtains the GP theory of dipolar BECs, 
which was exactly solved in the Thomas-Fermi (TF) ap- 
proximation [ll], [l2|. In order to go beyond the GP 
theory, we insert the Bogoliubov prescription into the 
grand-canonical Hamiltonian K = H — [j,N, with the 
chemical potential fi and the number operator N, and re- 
tain terms up to second order in the fluctuations. Then 
we decompose the fluctuating field operator according 
to <5-0(x) = YjJ [W„(x)a, y +V*(x)at]. Here, W„(x) and 
V v (x) are the Bogoliubov amplitudes, whereas dj, and 
a v denote bosonic creation and annihilation operators. 
The index v is associated with excited energy levels and, 
therefore, the sum excludes the v = mode, which rep- 
resents the condensate. By requiring the Hamiltonian to 
be diagonal, we derive the BdG equations 



L 3 jH UM (x,y) if w , v (x,y)W(y) 
J a V»,v(xj) %(x,y)Jiv„(y) 



e ( U ^) (3) 



where the diagonal matrix elements are given according 
to HuM^v) = <5(x-y)#Fi(y)+ **(yMnt(x-y)tf(x) 
and the off-diagonal ones according to Hjjy ( x , y) = 
^ r (y)Vi n t ( x — y) ^(x). In these expressions, an impor- 
tant role is played by the functionally diagonal part, 
which is given by the fluctuation Hamiltonian 

flpi(x) = h -fi + J dV**(x')M nt (x-x')*(x'). (4) 

Moreover, in the framework of the BdG theory, the oper- 
ator c\ v annihilates the ground state |0), so that excita- 
tions are interpreted as quasi-particles with energy e v . 

In general, solving the BdG equations exactly is very 
difficult and numerically arduous in the case of long- 
range interactions [24j]. Nonetheless, they provide an 
adequate starting point for investigating quantum fluc- 
tuations in most cases of interest. Typically, the number 
of particles is large enough, so that we can neglect the 
condensate kinetic energy and its wave function becomes 
vE'(x) = \/rJo(x), where rio(x) denotes the TF conden- 
sate density. Moreover, in this regime, a semiclassical 
approximation for the excitations is appropriate. Let 
q v (x) be either U v (x) or V v (x). The approximation is 
implemented through e v — > e (x, k), J2 U ~ * J d 3 fc/(27r) 3 , 
and q v — > q (x, k) e lk x , where g(x, k) is a slowly vary- 
ing function of x (25[. With this, the free Hamiltonian 
becomes ho(x, k) = h 2 k 2 /2M + C7tr(x). Despite the 
simplification, we still cannot solve Eqs. ([3]) because of 
their non-local character. The semiclassical approxima- 
tion also offers a procedure to deal with this problem. 



Consider the exchange term in Eq. ([3]), which, in gen- 
eral, reads /ex = / d 3 yHny (x, y) q v (y). Within the 
local density approximation (LDA), it becomes Iex ~ 
g(x,k)n (x)14 nt (k), with V- nt (k) = .g[l + e d d(3cos 2 0-1)} 
and cos 9 — k • e z . This is the leading order in a sys- 
tematic gradient expansion [261 ]. To estimate the error 
of ignoring the next order, we replace the gradients as 
V x — > 1/-Rtfj — > 1/K C , with i?TF the mean Thomas- 
Fermi radius and K c the momentum scale defined by 
the speed of sound. In terms of the condensate den- 
sity at the trap center, one finds the LDA to be valid 
for 3 x LY 2 a 3 n (0)]s x [1 + e dd (3 cos 2 9 - > 1. For 
edd < 1) this is commonly met in BEC experiments. 

One of the important physical features brought about 
by quantum fluctuations is the condensate depletion, i.e., 
particles are expelled from the condensate to excited 
states by interactions. The depletion density is, in gen- 
eral, given by n(x) — n (x) = Y^ v ^( x )V 1/ (x) and reduces 
in the present semiclassical theory to 



An(x) = ^Q 3 ( edd ; 



n(x)a 3 



(5) 



Here, n(x) denotes the total density of the gas and the 
auxiliary function Qi(x) = f} du (l — x + 3 x it 2 )^ 2 de- 
scribes the contribution of the DDI. In particular, the 
function Q3 is an increasing monotonic function of the 
relative interaction strength edd, which is bounded be- 
tween Q 3 (0) = 1 and Q 3 (l) = 3^/4 « 1.30. For 52 Cr, 
the most extensively studied dipolar system, the mag- 
netic moment of 6 /ib with the Bohr magneton /ib and the 
s-wave scattering length a s « 100 ao with the Bohr ra- 
dius ao yield the relative interaction strength e dd m 0.16. 
Then, the gas parameter n(x)a 3 evaluated at the trap 
center is typically of the order n(0)a 3 w 10~ 4 Q. For a 
homogeneous system with this value of the gas parame- 
ter, the fractional depiction amounts to AN/N « 1.5% 
and is not significantly affected by the DDI. 

The presence of quantum fluctuations also gives rise to 
a shift of the ground-state energy. The present BdG the- 
ory yields the following correction to the energy density 



AE(x) 
n(x) 



— gn(x)Q 5 (edd)y^^: 



(6) 



where Q5 is also a monotonic function and grows from 
Q s (0) = 1 up to Q 5 (l) = 3V3/2 w 2.60. Thus, phys- 
ical quantities directly related to the ground-state en- 
ergy offer much better perspectives for observing dipolar 
quantum fluctuations. Indeed, with the help of Eq. ((6|) a 
corresponding LHY equation of state for a homogeneous 
system can be derived, which leads to the Beliaev sound 
velocity 27 1 both in the absence [lU and in the presence 
of the DDI. Note that both functions Q3 and Q5 become 
imaginary for edd > 1- This arises from the fact that a 
homogeneous dipolar gas is only stable for edd < 1, oth- 
erwise the DDI dominates over the contact interaction 
and its attractive part leads to collapse. 
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Superfluid Hydrodynamics — At T = 0, the BdG the- 
ory states irrespective of the two-particle interaction that 
the whole sample is superfluid though not all particles are 
condensed. Thus, the system can be studied by means of 
the superfluid hydrodynamics. Indeed, by writing the su- 
perfluid velocity as v = V%, the superfluid hydrodynamic 
equations can be obtained by extremizing the action [28j 



A[n,x] 



dtd 6 xn I M 



+ en 



(7) 



with respect to the density n(x, t) and the velocity po- 
tential x(x, t). The energy density of the system as a 
functional of the particle density e[n] contains a mean- 
field component, which reads in TF approximation 

e M F = f/t r (x) + |n(x,t)+ /dV ^^ n^t), (8) 

and a quantum correction, which is given in the case 
of a dipolar Bose gas by the rhs of Eq. (j6]). Let us, 
now, extremize the resulting action ([7]) within a varia- 
tional approach. To this end, we adopt a harmonic ansatz 
for the velocity potential x(x, t) = [a x (t)x 2 + a y (t)y 2 + 
a z (t)z 2 ]/2 and a parabolic ansatz for the particle density 



n(x, t) 



15JV 

87Tif(i) 



f - 



V 1 



(9) 



for n(x,t) > and n(x, t) = elsewhere. The bar de- 
notes geometric average. Note that the particular density 
profile in (|9|) corresponds to the exact mean-field TF so- 
lution ll|, [l2j . This is consistent with the BdG theory, 
which assumes that interactions are not too strong. 

In order to derive quantum corrected equations of mo- 
tion for the TF radii Ri, we eliminate, at first, the varia- 
tional parameters for the velocity potential through their 
equations of motion on = Ri/Ri. Then, in the case of a 
cylinder symmetric trap, we obtain 



R.< 



-uj x R x - 



-lo 2 R 



15gN/4ir 



MR X R 



MR Z R 



I + 2e dd B 
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Rj' 
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~ 3 

R 2 



together with the auxiliary functions 
A(x) 



3 x 2 f s (x) . . 3 f s (x) 



2 x 2 - 1 ' v ~' 2x 2 — I' 

Here, f s (x) denotes the anisotropy function 

l + 2x 2 3x 2 tanh" 1 y/l - x 2 



fs(x) = 



1-x 2 



(1 - .T 2 ) 3 / 2 



(10) 



(11) 



(12) 



which arises for both BECs [111 |29| and nonsupcr- 
fluid fcrmionic dipolar gases [30(. In Eqs. (JTDJ), the 




c dd = 0, a s = lOOoo 
e dd = 0.16, a s = 100a 



Monopole 
Quadrupole 



A 



FIG. 1: (color online.) Quantum correction to the oscillation 
frequencies. 



quantum correction is governed by the parameter /3 = 
jQbi^d^^/a^N, with the constant 7 = 32 • 5^ • 7/2~ w 
4.49. Consequently, setting (3 = in Eqs. dTUll leads to 
the exact mean-field TF results 



ll|,ll2|. 



The beyond mean-field equations of motion for dipolar 
Bose gases fTO]) are the main result of the present letter. 
They allow for investigating both the static and the dy- 
namic properties of these systems. To this end, we solve 
them by considering /3 as a small perturbation, as this is 
consistent with the assumption of small fluctuations. 

Results - As is well known, the DDI affects the sys- 
tem in an anisotropic manner. Thus, in contrast to the 
case of a Bose gas with contact interaction alone, quan- 
tum fluctuations lead to a correction in the gas aspect 
ratio in equilibrium. To investigate this effect quanti- 
tatively, we set the rhs of Eqs. (fTTj|) to zero and solve 
them up to first order in j3. Thereby, one finds that the 
TF radii are given according to Ri = R® + 5Ri, where 
R° denotes the mean-field value and 6Ri represents the 
quantum correction. Thus, the beyond mean-field as- 
pect ratio is k = (R a x + SR X )/{R° Z + 8R Z ) w k° (1 + 6k). 
Usually, the correction to the aspect ration in the present 
experiments is quite small. For 52 Cr for example, one 
has at most 8k < 1 0~ 2 , which is too small to be detected. 
For stronger interactions, this turns out to be an impor- 
tant effect, since it implies a characteristic dependence 
on the trap geometry due to the DDI. Consider, e.g., the 
experimental values for particle number N — 3 x I0 4 
and trap frequencies Zu ps 2tt x 660 Hz as in Ref. 0], 
but instead of chromium values, use a s = 190 ao and 
£dd = 0.7. This could represent a dysprosium sample, 
which is heavier than chromium and possesses a larger 
magnetic moment of 10 /Kb yielding the dipolar length 
add ~ 133 ao- Then, the correction for the aspect ratio 
would be about 5k w 27 % for moderate cigar-like traps. 

Let us, now, analyze the low-lying excitations. To 
this end, we apply for the TF radii the ansatz Ri(t) = 
Ri(0)+r)i sin (ftt + <f) and obtain by linearizing Eqs. (fT0|) 
a matrix equation which leads to the frequencies of the 
monopole and of the quadrupole modes up to first order 
in /?. In Fig. HJ we plot the resulting corrections to the os- 
cillation frequencies as functions of the trap aspect ratio 
A in units of 50. = (63^/128)^03 n(0) for e d d = 0.16. 
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FIG. 2: (color online.) Aspect ratio as a function of time. 

The dashed lines represent the Pitaevskii-Stringari result 
for contact interaction alone [HI , which we immediately 
recover by setting edd = 0. Fig. [1] shows a dependence 
on the trap aspect ratio, which is substantially different 
from the contact case and stems from the quantum cor- 
rection for the aspect ratio. Thus, tuning the trap aspect 
ratio could be used to identify the effects of the DDI 
at the many-body level in ultracold Bose gases. This, 
however, requires stronger interactions. Indeed, the os- 
cillation frequency of the intermediate mode of 52 Cr has 
recently been investigated in a triaxial trap [H . For that 
experiment we estimate the quantum correction of the 
oscillation frequencies to be of the order SCI « 3 x 10~ 3 , 
due to the small particle number and weak trap. For 
more favorable values of these parameters, such as those 
in Ref. Q, one already obtains <5X2 ps 1 %. For the previ- 
ous hypothetical dysprosium setup, one could even have 
SCI w 4 %, which should render these effects observable. 

For the TOF dynamics, the quantum correction of the 
equilibrium TF radii also plays an important role. In- 
deed, by inserting the expansion Ri(t) = Ri(t) + SRi(t), 
one derives a set of coupled differential equations for the 
mean- field TF radii Ri(t) and their corrections SRi(t). 
For a system consisting of 52 Cr atoms, no signal of quan- 
tum fluctuations can be seen in the present TOF analysis 
and the mean-field theory provides again a good agree- 
ment with the experimental results For stronger in- 
teractions, however, the situation changes. In Fig. [21 we 
plot the ratio «(t) = [R° x (t) + 8R x (t)]/[R°(t) + SR z (t)} 
as a function of time for the previous dysprosium setup. 
Then, clear differences with respect to the mean-field 
TOF dynamics show up, specially at large times. To 
investigate the dependence of the asymptotic aspect ra- 
tio on the system parameters we have settled the dipo- 
lar length at add ~ 133 ao and have varied the s-wave 
scattering length o s , so as to vary the relative strength 
from edd = 0.2 up to edd = 1-0. Our results are 
shown in the inset of Fig. [21 where the asymptotic value 
[k(oo) — k°(oo)]/k°(oo) is plotted against edd for differ- 
ent trap aspect ratios. The results reasonably deviate 
from the mean-field values, so that quantum fluctuations 



should be observable in the TOF dynamics of Bose gases 
of highly magnetic atoms. 

Conclusion - Applying the BdG theory in the limit 
of large particle numbers, we obtained analytic expres- 
sions for the condensate depletion and the ground-state 
energy of dipolar Bose gases. This allowed us to study 
key properties of the system such as low-lying excita- 
tions and TOF dynamics beyond the mean-field approx- 
imation. We find that both polar molecules and highly 
magnetic atoms are realistic scenarios for observing un- 
precedented beyond mean-field dipolar physics. 

We thank J. Dietcl and L. Santos for useful discussions 
and the DFG for financial support (KL256/53-1). 
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